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Abstract 

We perform ab initio QED calculations of energy levels for the n = I and n = 2 states of He-like 
ions with the nuclear charge in the range Z = 12-100. The complete set of two-electron QED corrections 
is evaluated to all orders in the parameter aZ. Uncalculated contributions to energy levels come through 
orders a^(aZ)^, a^{aZy, and higher. The calculation presented is the first treatment for excited states 
of He-like ions complete through order a'^{aZ)^. A significant improvement in accuracy of theoretical 
predictions is achieved, especially in the high-Z region. 

PACS numbers: 12.20.Ds, 31.30.Jv, 31.10.+Z 
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Introduction 



Helium and helium-like ions, being the simplest many-electron systems, traditionally serve as 
an important testing ground for investigations of many -body relativistic and QED effects. Cal- 
culations of QED effects in He-like ions have a long history. The expression for the Lamb shift 
complete through orders a{aZy and a'^(aZ)^ was derived in pioneering studies by Araki [ 1] and 
Sucher [Q]. Numerous posterior investigations of higher-order QED corrections in two-electron 
systems (see, e.g., review |||3|fl and recent original studies fllsiyj) were primarily aimed at helium, 
in which the experimental accuracy is by far better than in other two-electron systems. Recent 
progress in experimental spectroscopy of highly charged ions |Q. S 1^ opened new perspectives 
for probing higher-order QED effects in ions along the helium isoelectronic sequence up to He- 
like uranium. Investigations of QED effects in high-Z ions are of particular importance since they 
can provide tests of quantum electrodynamics in the region of a very strong Coulomb field of the 
nucleus. Another factor that stimulates these investigations is the possibility to test the standard 
model by studying the effects of parity non-conservation (PNC) lllflL Experimental 
identification of the PNC effects will require precise knowledge of the 2^5*0 — 2^Po interval in He- 
like ions with nuclear charge numbers near Z = 6A (gadolinium) and Z = 90 (thorium), which 
happens to be very small for these values of Z thus enhancing the PNC effects significantly. 

Investigations of QED effects in heavy He-like ions differ significantly from those for the 
helium atom. First of all, the nuclear coupling parameter aZ approaches unity and cannot be 
regarded as a good expansion parameter as in the case of helium. But on the other side, the 
electron-electron interaction in these systems is suppressed by a factor of 1/Z with respect to the 
electron-nucleus interaction and, therefore, can be accounted for by a perturbation expansion in 
the parameter 1/Z. 

Until recently, the only QED effects calculated to all orders in aZ were the one-electron self- 
energy and vacuum-polarization corrections ll^ [isl . So, theoretical investigations of energy 
levels in heavy He-like ions mostly relied on these one-electron values, correcting them to account 
for the "screening" effect by various semi-empirical rules, notably, within Welton's approximation, 
as in Ref. Q. A more elaborate treatment of QED effects in He-like ions was presented by Drake 
Iryj . His values for the QED correction included the complete contribution to order a'^{aZ)^ 



derived in Refs. Ill, 121] and parts of higher-order contributions obtained by employing the all-order 
results available for the one-electron QED corrections. The total energy values of Ref. 1 17] are 
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complete through order a'^(aZY and uncalculated terms start in orders a^(aZ)^ and a^{aZY. 

Later, Johnson and Sapirstein lisl applied relativistic many-body perturbation theory (MBPT) 
to the treatment of the electron correlation forn = 2 triplet states of He-like ions. Combined with 
Drake's values for the QED and recoil corrections, their results yielded a better agreement with the 
experimental data than those of Ref. [ 17]. While the approach of RefJ 18|] is still incomplete to 
order a^(aZ)^, it includes terms that were not accounted for in Ref. flTj, namely the Breit-Breit 
interaction and some relativistic corrections to the second-order energy. Later, other evaluations 
of the electron-correlation part of the energies of He-like ions were performed by the relativistic 
configuration-interaction (CI) method lliol] and by the relativistic all-order MBPT approach OOl l. 

nnri 

The studies II18L I19L 12011 share the same main features: their treatment is based on the no-pair 
Hamiltonian and the electron correlation is taken into account within the Breit approximation. 
The results of these evaluations are in a very good agreement with each other. 

A somewhat different approach was employed in Refs. 112 iL 12211. While the electron-correlation 



part was evaluated (as in the prev ious work by the same group UM) by the CI method, the QED 
part was not taken from Ref. |17f| but evaluated independently, by considering the one-loop QED 
corrections in a local screening potential. Due to different treatments of QED effects, there are 



certain deviations between the results of Refs. ||21 



2M and those of Refs. ||18l ll9ll2 



mplished 



20,l2i|]- To perform similar 



In order to obtain reliable predictions for energy levels of high-Z ions and to improve the 
theoretical accuracy in the low- and middle-Z region, it is necessary to take into account two- 
electron QED effects without an expansion in aZ. Such project has been recently accornplished 

ai to order a^) for the two-electron part of the ground-state energy of^ He-like ions 
26l 1271] and for the lowest-lying states of Li-like ions 11281 
QED calculations for excited states of He-like ions is more difficult. One of the reasons is that, 
for the first time in QED calculations to all orders in aZ, we encounter levels that are quasi- 
degenerate, namely 2^ Pi and 2^ Pi. To derive formal expressions for QED corrections in case of 
quasi-degenerate states is a serious problem that has been solved first within the two-time Green 
function (TTGF) method 11321 33L 3411 . Different approaches to this problem have recently been 
addressed by other authors psLbq] . 

Several QED corrections have been calculated to all orders in aZ for excited states of He-like 
ions up to now. In our previous investigation Irk, we evaluated the vacuum-polarization screening 
correction for all n = 2 states of He-like ions. The two-photon exchange correction was calculated 
for excited states of He-like ions by Mohr and Sapirstein i38ll (2^ Si and 2^Po,2 states), by Andreev 
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(2^'3Pi), and by Asen et al. (2^S'o,i). In this paper we 
present an evaluation of the self-energy screening correction and an independent calculation of the 
two-photon exchange correction for all n = 2 states of He-like ions. This completes the ab initio 
treatment of all two-electron QED corrections of order to all orders in aZ and significantly 
improves the theoretical accuracy for the energy values, especially in the high-Z region. Unlike 
all previous calculations, the results obtained are complete through order a^{aZY; uncalculated 
terms enter through three-photon QED effects (to order a^{aZY and higher) and through two-loop 
one-electron QED corrections {a^{aZy and higher). 

The paper is organized as follows. In the next section we describe the basic formalism and 
present general formulas for the two-electron QED corrections for the case of quasi-degenerate 
levels. In Section lUl the numerical procedure is briefly discussed and numerical results are pre- 
sented for the two-photon exchange correction and the screened self-energy correction. The total 
two-electron QED correction is then compiled, analyzed, and compared with the known terms of 
the aZ expansion. In the last section, we present a compilation of all contributions available to the 
energy levels and compare results of different theoretical evaluations with existing experimental 
data. The relativistic units (/i = c = m = 1) are used throughout the paper. 



et al. 



(21^0 1, 23Po) and 136 



I. FORMAL EXPRESSIONS 



A. Basic formalism 



In this section we briefly formulate the basic equations of the TTGF method for quasi- 
degenerate states of a He-like ion. A detailed description of the method and, particularly, its 



ana, partic 

^. fly, 



m- The 



implementation for the case of quasidegenerate states can be found in Refs. 
derivation will be given for two particular quasidegenerate states, (ls2pi/2)i and (Is2p3/2)i, and 
can immediately be extended to a more general case. The unperturbed two-electron wave functions 
in the jj coupling are given by 

ui= ^ Oam,j,m,|JM)-^5^(-l)^|PaPt;), (1) 

mairiv ^ P 

U2= ^ Oam„j^m^|JM)-^5^(-l)^|PaPw), (2) 
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where a, v and w are taken to represent Is, 2pi/2 and 2^3/2 orbitals, respectively; P is the permu- 
tation operator: 

^(-l)^lPaPt;) = \av) - \va) , 
p 

and \av) = \a)\v) is the product of the one-electron Dirac wave functions. The transition to the 
wave functions corresponding to the LS-coupling scheme within the non-relativistic approxima- 
tion can be performed by 



with 

1 f V2 -1 , 

We mention that this choice of the matrix R implies that the one-electron 2pi/2 and 2^3/2 wave 
functions have the same sign in the non-relativistic limit. 

The standard definition of the four-time two-electron Green function in the external field of the 
nucleus is 

, X2', Xi, X2 ) = (0|TV^(x;)V^(x'2)^(xi)V^(x2)|0) , (5) 

where iJj{x) is the electron-positron field operator in the Heisenberg representation, ijj = ^^^7°, and 
T denotes the time-ordered product operator. This Green function is constructed by perturbation 



theory after the transition to the interaction representation where it is given by (see, e.g., 11431] ) 



, , X (0|T?/'in(x'i)^in(a;'2)^in(a;2)^/'i„(a;i)exp [-i/d^2;Hint(2:)] |0) 
GlXi , Xn: Xi, X2) = — F ? — : , . . (6) 

Here ipm{x) is the electron -positron field operator in the interaction representation and Hint is the 
interaction Hamiltonian. Expression Q allows one to construct G by using Wick's theorem. 

In what follows, it is more convenient to work with the Green function in the mixed energy- 
coordinate representation, which is defined by 

G'(pf,x;,p;°,x'2;p?,xi,p°,X2) = dx[dx'2dx[Ux'2" 

X exp {ipfxf + 1^2X2 — ip\x\ — ip^^^) G{x\, x'2] Xi, X2) ■ 

(7) 
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The Feynman rules for G{pi, x'^ , Xg; p?, xi , p^, X2) can be found in js^El- We now introduce 
the Green function g{E) as 

poo 

g{E) 6{E -E') = - dpi dpi dp[' dp',' 5{E - - p°) 

x5{E' -p'^ -p'^)PoG{ptp'^:plpl)lhlPo, (8) 

where Pq = '^ku\ projector on the subspace of the unperturbed quasi-degenerate states 
under consideration [see Eqs. dU) and Q]. It can easily be shown (see, e.g., Refs. QB) that 
the Green function g is the Fourier transform of the two-time Green function projected on the 
subspace of the unperturbed quasi-degenerate states. 

It can be derived (see Ref. |34] for details) that the system under consideration can be described 
by a two-dimensional Schrodinger-like equation (A; = 1, 2), 

H^k = Ekipk , ^k' = hk' , (9) 

where 

H = p-^/^KP-^/^, (10) 

K = ^^j^dEEg{E), (11) 

P = ±-J^dEgiE), (12) 

r is a contour in the complex E plane that surrounds the levels under consideration but does not 
encircles other levels, and E^ are the exact energies of the states under consideration. It is assumed 
that the contour T is oriented anticlockwise. The operator H, which is a 2 x 2 matrix, is constructed 
by perturbation theory in a. Substituting 

g{E) = g(^)(E)+g('\E)+g('\E) + ..- , (13) 
P = p(0) + pW + p(2) + . . . , (14) 

K = + K^'^ + K^'^ + ■■■ , (15) 



where the superscript indicates the order in a, we obtain 113 3h 



= , (16) 
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^ ^(2) _ 1 p{2)^(0) _ 1^(0) p(2) _ 1 p{l)^(l) _ 1^(1) p(l) 

2 2 2 2 

+^ p(i) p(i) /^-w I ^7^(0) p(i) p{i) , 1 pWj^m p(i) /I ON 

8 8 4 ■ ^ ^ 

The solvability of Eq. ^ yields the basic equation for the calculation of the energy levels 

det{E~H) = 0. (19) 



As was noticed in Ref. 113 311 . due to nonzero decay rates of excited states, the self-adjoint part of 



H should be understood in Eqs. ^ and ( [T9t . 

H={l/2){H + H^). (20) 
To zeroth order in a, the Green function g{E) is 



2 



s=l ^ 

where Ef^ and Ef^ are the unperturbed energies of the (ls2]9i/2)i and (Is2p3/2)i states, respec- 
tively, given by the sum of the one-electron Dirac-Coulomb energies: 

i-,(o) _ I ,-,(0) _ I 

-C/l - eis+ £^2pi/2 ' £/2 - £ls + £^2p3/2 • 

Substituting Eq. (|2T1) into the definitions of K, P, and H, one gets 

K^=Er^5^,, (22) 
4°^ = (23) 
=i^r^..- (24) 

Now we introduce a set of notations that will shorten the following expressions. The short-hand 
notation will be used for the summation over the Clebsh-Gordan coefficients in Eqs. ([T]), Q: 

Fi\iii2) = {jhmiJi2mi2\JM)\iii2) . (25) 

where \i1i2) is either \av) or \aw). It is convenient also to use the notation for the operator of the 
electron-electron interaction: 

JH = e2<«2'^/..H, (26) 
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where = 7^7^ = (1, o:) and D^^, denotes the photon propagator. In the Feynman gauge, the 
propagator of a photon with the non-zero mass fi is 



exp ijJ^ — 


/i2 + zO 


|x-y|) 


47r 


|x- y| 





B^,^{u), X - y) = g^^ — ■ , (27) 



where it is assumed that \m.^ u"^ — + zO > 0. For the matrix elements of the operator we 
will use the short-hand notation 

h.uiyj) = ■ (28) 



B. One-photon exchange diagram 

In order to illustrate how the method works, below we present the detailed derivation of the 
correction to the quasidegenerate energy levels (ls2pi/2)i and (1 52^3/2)1 due to the one-photon 
exchange diagram (Fig. [l]). While the corresponding evaluation is much less cumbersome than 
those for the second-order two-electron corrections, it demonstrates most essential features that 
are encountered in these cases. For simplicity, in the derivation below we will assume that the 
unperturbed energy of the initial state i differs from that of the final state k: Ef^ 7^ Ef^ (in the 
case under consideration it corresponds to z 7^ k). However, all the final formulas can be shown to 
be valid also for the case Ef^ = Ef^ . 

According to the Feynman rules and the definition of g{E), the contribution of the 

one-photon exchange diagram is 

^ IpiiPi2kik2{p'i -Pi) .29) 

{p^,-ek,+tO){E-p^,-Sk,+iO)' 

Employing the identities 

1 1/1 1 



{pf - epi, + 20) {E - pf - Bp,, + zO) E- e'^^ \Pi - ^Pi. + ^0 E- pf - + zO , 

(30) 



+ n ^ , (31) 



(p? - + i^){E -p\- Sk, + ^0) E- Vrf - ^fci + «0 E-p\- ek, + , 
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we obtain 



(1) 

ik 



FiFk— i dE 



E 



27ri/r {E-Ef^){E-El 
1 1 



X 



p'l" - epn + ^0 E- - epi, + iOj VPi -Sk.+tO E-p\- e^, + iO 



+ 



1 



'xIphPi2k1k2i.P1 — Pi) 



(32) 



The expression in the square brackets is an analytical function of E inside the contour F, if the 
photon mass /i is chosen properly (see Refs. Carrying out the E integration by Cauchy's 

theorem and taking into account that 



1 1 

+ 



2tt I \ x + zO —X + zO 



5{x) 



we obtain 



pEf^ IphPi2kik2{^Pii - Pi) 



^(0) _ ^(0) 



X 



+ 



I 



p\-ek,+i^ E^^ - p1 - ek2 + iO 

/oo 
dp?J2{-i) 
-00 n 



pEfhpi^Pi^kik2{pi - £ki] 



^(0) _ ^(0) 

1 



X 



p'r-ep,,+tO E^^^ -pf-ep,,+iQ 



In the same way we find 



.pIphPi2kxk2{^Pii ~ Pi) 



+ 



p\-eu,+tQ Ef^ -p\-ek2+tQ, 

1 / , /n X ^/ ^ \pIpi\Pi2k1k2 [P'l - ^k-i) 



+ 



(0) _ p(0) 

1 



Er-E 



pf - epi, + ^0 4°) _ p'o _ ^^.^ + ,0 



(33) 



(34) 



(35) 



9 



Substituting Eqs. (l34b . (l35t into Eq. (flTt . we get 



X — r + 



■ POO ^ 



P 

1 



- ^Pn + ^0 + _ ^,0 _ + ,0 j J ■ ^^^^ 
Introducing the notations Ai = £pi^— and A2 = Spi.^— Sk^, we can rewrite Eq. (B^ as follows, 

/OO / -I -1 

-00 p + ^1 + *0 A2 - + 20 

1111 
+ 1 + + 1 + 



+ A2 + Ai - + iO - Ai + iO -A2 - uo + iQ 

1 1 

+ 



(37) 



u;-A2 + iO -Ai-cu + iO 

= ^i^k^ / y~l(~l)^-^PnPi2fciA:2('^) 

^ J —CO p 

X + Ai) + - Ai) + 5{uj + A2) + 5{uj - A2 
Taking into account that I{z) = I{—z), we finally obtain 

H-l^ = PiPk- — l)'^[-^PiiPj2A:iA;2(^l) + -^PjiPi2fciA;2 (^2)] • (38) 
P 

C. Two-photon exchange diagrams 

The set of two-photon exchange diagrams is shown in Fig. El The first and the second graph 
are referred to as the ladder and the crossed diagram, respectively. The derivation of the general 
expressions for the two-photon exchange correction in the case of quasi-degenerate levels is rather 
lengthy. However, it greatly resembles the corresponding derivation for the one-photon exchange 
correction presented above, on one hand, and that for the two-photon exchange diagram in case 



of a single level described in detail in Ref. ll45h . on the other hand. We thus present only the final 



formulas for the two-photon exchange contributions to the matrix elements of the operator "^^^ . 
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1. The ladder diagram 



The contribution of the two-photon ladder diagram is conveniently divided into the irreducible 
and the reducible part. The reducible contribution is defined as a part in which the total inter- 
mediate energy of the atom equals to Ef^ or E^2^ and the irreducible part is the remainder. The 
operator H^"^^ is defined by Eq. dTSt . The first three terms in the right-hand side of this equation 
contribute both to the irreducible and to the reducible part. As to the others, it is natural to ascribe 
them to the reducible part. 

The contribution of the irreducible part of is defined as the self-adjoint part of the follow- 
ing matrix 



The result is 

\s,k{Ef\ 0, 0) + S,k{Ef\ 0, A) + 0, 0) + S:,k{Ef\ -A, 0) 



1 

4 L 



? r°° 1 r 
+ — v.p./ dx- S,k{Ei'\0,x)-S,k{Ei'\0,x + A) 
47r J^^ X L 



+S,k{E^^\x,0)-S,k{E^^\x-A,0) 



where A = Ef^ — E^^^ and the matrix elements Sik are defined by 



(40) 



SikiE,x,y) = ^{-1 



2tt 



X 



°° , Sr^ IphPi2nin2{^Pii — ^ + X) /nin2fcifc2(^fci —^ + 11) 



[UJ - Sn^ (1 - ZO)] [E -UJ - (1 - 



The summation here runs over all ni and for which E'i^ ^ Ef\Ef\ where = + e. 
is the total intermediate energy of the atom. The sign "v.p." in front of the integral in Eq. (HUl) 
denotes that the principal value of the integral (over x) must be taken. 

We note that the part containing the integral over x in Eq. (l40t vanishes identically in case 
of diagonal matrix elements {i = k). It neither appears for single levels 1I45I1 . In case of off- 
diagonal matrix elements (i 7^ k), the contribution of this part is of order a^A and it vanishes 
when {Ef^ - E^) ^ 0. As shown in Ref. [Q], such terms contribute to the next order of 
perturbation theory and can, therefore, be disregarded in the present consideration. Expression 



n2 
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(l40b can be simplified even further by taking into account that 
where = {Ef^ + Ef'^)/2. We thus write '^^ simply as 



lad, ir 
ik 



• /.oo 



p(0) , p(0) p{0) 



X 



EI Pi\Pi2n\n2i,^ Pi\ "-i^) -^nin2feifc2 (^fci ^) , ^/ 2 an 
~ — 3;;;;=^ ' ^ + ^^"^ ^) • 



nma [cu - - «0)][E - - (1 - iO)] 
The reducible contribution is induced by the self-adjoint part of the following operator 



lad, red ^lad,rcd,a _j_ ^-lad, rcd,6 



where 



and 



^lad,red,a _ ^(2,red) _ ^p(2, red)^(O) _ ^^(0) p(2, red) 

2 2 



^lad,red,6 ^ _lp{l);^(l) _ l;^(l)p(l) , ^p(l)p{l)^(0) 

2 2 8 

I ^7^(0) p(l) p(i) , 1 p(l)l^(0) p(l) 
8 4 



The result for the first part reads 



H 



lad, red, a 
ik 



ikl 



D,k{Ef\ 0, 0) + Afe(^r, 0, A) + Afe«^ 0, 0) + Afc(£^r, 

(0) n n..^ 77^(0) 



.(0) 



.(0) 



.(0) 



--^v.p./" dx- Afc(^r,o,x)- Afc(^r,o,x+A) 



+D,k{E^^\x, 0) - D,k{El^>,x - A, 0) 



(0) 



where 



p(0)_p(0) p(0) 



I Pi\Pi2n\n2 



^ni ) ■Inin2k 



ik2 i^ki 



-ni > 
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Bik{x) 



-fin — -t'l i-t'2 



X 



oo (lU 



I Pi\Pi2n\n2ij^ Pi\ ^ni ) -^711712^1^2 ^ 



(0) 



zO)(to' — efci + X — iO) 



(49) 



p(0)_p(0) p(0) 



p 



27r 



nin2 



Tlx, 



—00 
00 



- .0) 



27r - - z0)(a7 - + Ef - - zO) 



(50) 



r(0)_p(0) p(0) 



27r 

P n,in2 

r-oo Ipi^Pi^riYn^k^Pix — + ^c) In^n^k^k^K^kx ^ ^ + U) 

X / da; 7 -7 r 

[UJ - - %^){bJ + en2- E - tO) 



(51) 



The part containing the integral over x in Eq. (l47l represents a contribution of order A. Again, 
we regard this contribution as belonging to the next order of perturbation theory and disregard it 
in the present investigation. 

The second part of the reducible contribution is given by the matrix element of the operator 
(l46b . The result is obtained by taking into account that 



^ik - Ei Oik , 



p(0) _ r 
Uk - Oik , 



(52) 



K. 



(1) 

ik 



FiFk E](~l)^ I 2 [■^PiiPi2kik2{.^l) + IpiiPi2kik2{.^1 

{Ef+ET) ^ 



2Ti 



du I pi-^Pi^kxk2{^) 



_(a; + Ai - iQ){uj - A2 - iO) (u; + A2 - iQ){uj - Ai - 'iO)_ 



(53) 



and 



P. 



(1) 



ik 



-FiFk / (—1) — / du Ipi^pi2kik2i^) 



+ 



a; + Ai - iO){uj - A2 - iO) (u; + A2 - iO)(a; - Ai - zO)_ 



(54) 



13 



The total result for the reducible part can be simplified by using Eq. (l42l) and disregarding terms 
that contribute to the next order of perturbation theory. One can show that in this case the A's, 
B's, and C's in Eq. (H71) are cancelled completely by the ijiad,rcd,b jgj^ jj^g result is just 



p{0)_p(0) p(0) 



2tt 



P nxn2 

-?0)(c^ + £„, -E^ ^ -zO) 



2. The crossed diagram 



The contribution of the crossed diagram is induced by the self-adjoint part of the following 
operator 

^cr _ ^(2) _ (l/2)p(2)i^{0) _ (1/2)K(°)P(2) _ (56) 

The corresponding result reads 



rrcr 









i 






+— 


v.p. 


[ dx — 


47r 


'-00 a; L 



.(0) 



.(0) 



n(0) 



T,fc(Ef \ 0, a;) - T,fc(Ef , 0, X + A) 
+Tik{E^^\ X, 0) - T,fc(Ef , X - A, 0) 



(57) 



where 



fOO 

X / (iu; ■ 

-oo 



nin2 



J- Piin2nik2 

[epi^ -UJ + x) 



■ (58) 



[uj - En^il - iO)][E - Epi^ -Sk^-x-y + uj - £:„2(1 - zO)] 

The expression ( I57l can be simplified in the same way as the previous contributions, with the 
result 



HTk = FiFkT,k{E^'\0,0) + O{a^A)=F,FkJ2(-^fYl 

P nin2 
Ipiin2nik2{^Pii '~ ^) -^niPi2fcin2 (^^i ~ 



2^ 



fOO 

X / duj 

-oo 



[U - Eniil - iO)][E - Epi^ - Ekj + - - iO)] 



+ 0{a^A) . (59) 
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D. Screened self-energy correction 



The set of Feynman diagrams representing the screened self-energy correction is shown in 
Fig. |3] Formal expressions for this correction in case of quasi-degenerate states were obtained 
previously in Ref. MoA by the TTGF method. Here we present only the final expressions for this 
correction. 

The contribution of the vertex diagrams is given by 



= F^F, 5:(-i)^f r du y: — ^^i^^ 



+ 



711712 

Piinikin2 

(A2) Ip 



iO)][ek^ - u - en2{l - iO)] 



(60) 



[epi2 - u - en^{l - i0)][ek2 - u; - £:n2(l - ^0)] 
where Ai = epi^ - Ek, and A2 = epi^ - ek2- 

The contribution of the remaining diagrams is conveniently separated into the irreducible and 
reducible parts. The irreducible contribution is given by 

'PilPj27lfc2 
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n^Pi2 
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where is the self-energy operator defined by its matrix elements, 



(a|S(e)|6) 



27r 



/oo 
-00 „ 
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The result for the reducible contribution reads 



Hir^ = F^Fk \ 5^(-l)^{/p.,«,fc,fc,(Ai) [(P^i|S'(e«J|Pzi) + (fci|S'(5,J|fci) 
p 

+IphPi2kik2{^2) {Pi2\^'i£pi2)\Pi2) + (A;2|S'(£:fcJ|/c2) 
+I'p.,P.,k,k2i^i)[{Pn\nep.J\Pt^) - {h\j:i6k,)\h)_ 
+I'pnPr2k,k2i^2) [{P^2\nep^2)\P^2) - {k2\nek2)\k2)\ ] + 0{a^A) 
where = dl{uj)/duj, and S'(cu) = dT.{u)/duj. 



(61) 



(62) 



(63) 
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E. Screened vacuum-polarization correction 



The derivation of formal expressions for the screened vacuum-polarization correction in case 
of quasi-degenerate states was described in our previous work IstI . For completeness, we present 
here the final expressions for this correction; the corresponding set of Feynman diagrams is shown 
in Fig. 121 

The expression for the contribution of the diagram with the vacuum-polarization loop inserted 
into the photon propagator can be obtained from the formula for the one-photon exchange (l38t by 
replacing the operator of the electron-electron interaction I{e) by the modified interaction, 



[/^p(£:,x,y) = ^ I duj I rfzi rfzg 
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xTr KG(c<; - e/2, zi, za) a" G{uj + e/2, zs, zi)] 



(64) 



where G{ijJ, x, y) = Ylm V^n(^)^ri(y)/[^ ~ ^n(l — ^0)] is the Dirac-Coulomb Green function. The 

(2) 

corresponding contribution to iJ^-^ is 



o-vp, ph 



PllPl2\U^p{^l)\hk2) + {PtlPl2\U^l{A2)\hk2) 



(65) 



where Ai = epi^ - Sk^ and A2 = Spi^ - Sk^- 

To the order under consideration, expressions for the remaining diagrams can be obtained from 
the one-photon exchange correction by perturbing the wave functions and the binding energies by 
an additional vacuum-polarization interaction. The result is 



TT-VPjWf 



1 



^^vp,be ^ F,Fk-J2{-'^f[{SP^lP^2\[I{Ai)+I{A2)]\hk2) 
P 

+ {P^^6P^2\[IiA,) + I{A2)]\k^k2) 

+ {PtlPt2\[I{Ai) + I{A2)]\6kik2) 
+ {Pt,Pt2\[I{A,) + I{A2)]\k,6k2) 

+{6ep,, - fe,J {Pt,Pt2\I'{A,)\k,k2) 
+{5epi, - fefcj {PiiPi2\I'iA2)\k^k2) 



where 5i and 5k refer to the first-order corrections to the corresponding wave function. 



(66) 



1^^) = E 



\n){n\\J\jTp\i) 



(67) 
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5ei is the correction to the energy, 5ei = {i\Uyp\i) , and 

t/vp(x) = — duj dy- ^Tr [G(^, y, y)] (68) 

^TTi j_oo J |x-y| 

is the vacuum-polarization potential. 

As discussed previously in Ref. a direct derivation based on the TTGF method yields a 
result that differs from Eq. (l66b by terms of order (a^A), which can be disregarded as long as we 
are not interested in higher orders of perturbation theory (see Ref. I34I] for a detailed discussion). 

II. NUMERICAL EVALUATION AND RESULTS 

An important difference of the present investigation from the previous studies of QED effects 
in high-Z ions is that it involves QED corrections for quasi-degenerate configurations, namely 
(ls2pi/2)i and (ls2p^/2)i. While the derivation of basic expressions in this case is more difficult 
than for a single state, the final expressions for the diagonal matrix elements turn out to be very 
similar to those for the single-level case. We can, therefore, adopt a code developed for single- 
level calculations for the diagonal matrix elements of the operator H. For an evaluation of the 
off-diagonal matrix elements, a generalization of the code is needed. 

The numerical procedure employed in the present calculation of the two-photon exchange cor- 



47|] 



rection is based on that presented in detail in our previous investigations for Li-like ions ||30 
Apart from the angular reduction that is performed by using the standard angular-momentum tech- 
nique, the evaluation is rather similar to that for Li-like ions. The calculation was carried out em- 
ploying the Fermi model for the nuclear-charge distribution, with the nuclear charge radii specified 
in SectionHni The numerical uncertainty of the results is expected to be 1 x 10^"^ eV in all cases 
except for the off-diagonal matrix element, for which the uncertainty is 1 x 10^^ eV for Z < 50, 
2 X 10~^ eV for Z < 80, and 4 x lO^'' eV otherwise. As a check of the numerical procedure, we 
performed the evaluation in two different gauges, the Feynman and the Coulomb ones. The two- 
photon exchange corrections (for mixing configurations, individual matrix elements) were found 
to be gauge invariant well within the uncertainty specified. 

The results of our numerical calculation of the two-photon exchange correction for n = 1 
and n = 2 states of He-like ions are presented in Table U The values listed represent correc- 
tions to the energy in case of single levels and contributions to the matrix elements for the 
quasi-degenerate states. The energy levels for the (ls2pi/2)i and (Is2p3/2)i states are obtained 
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by diagonalizing the 2x2 matrix H containing all relevant corrections. In Table HI we present 
also a comparison of our numerical values with the results of t he p revious calculations of this 



36 



38, 



or tnepi 

HQ. 



mm 

14 ID . The comparison indicates 



correction for various states of He-like ions II23L 
that calculations by different groups are generally in agreement with each other. However, there 
exist also certain deviations between different calculations, notably with those by Andreev et al. 

HQ. 

Regarding the comparison of the present results and the ones of Ref. Il36ll for the mixing 
states, we would like to stress that, generally speaking, results of different methods for individual 
matrix elements could be different, since the matrix H can differ by a unitary transformation. We 
observe, however, that in our case the results for the individual matrix elements agree with those 
of Ref. Isel approximately at the same level as for the single states. 

The calculation of the screened self-energy correction for n = 2 states of He-like ions resem- 
bles that for Li-like ions described in our previous work [29]. A more difficult angular structure 
of the initial- state wave functions for He-like ions makes final expressions more lengthy and their 
numerical evaluation more time consuming. Significant complications appear in performing angu- 
lar integrations in momentum space for the vertex part with free-electron propagators. To handle 
them, we developed a generalization of the angular- integration procedure described in Ref. W to 
arbitrary states, using our experience in calculating similar angular integrals for the two-loop self- 
energy diagrams MM- The actual calculation was carried out employing the spherical- shell model 
for the nuclear-charge distribution. Our numerical results for the screened self-energy correction 
for n = 1 and n = 2 states of He-like ions are presented in Table Hll in terms of the dimensionless 
function F(aZ) defined as 

AE = a^{aZfF{aZ) . (69) 

The values listed in the table represent corrections to the energy in case of single levels and con- 
tributions to the matrix elements Hik for the quasi-degenerate states. 

In case of the ground state of He-like ions, the self-energy correction was evaluated previously 

n n o 

by Persson et al. ||25I], by us |I27I1. and by Sunnergren ||49(]. In the present work, we recalculated 
this correction using the new code and found an excellent agreement with our previous results and 
with those by Sunnergren. A small deviation of the present result for Z = 90 from the old one is 
due to a more recent value for the nuclear charge radius used in this work. 

We note that the values presented in Tablellllfor n = 2 states of He-like ions can also be used for 
determining the screened self-energy correction due to the interaction of the valence electron and 
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the (Is)^ shell in Li-like ions. Indeed, by using elementary angular- summation rules, we obtain 

(2j. + 1) /\E];' = Y.i'^J + 1) AEf^j , (70) 
J 

where Ai?^^ denotes the screened self-energy correction in a Li-like ion due to the interaction 
of the electron in the state v and the (Is)^ shell, AE^^j is the screened self-energy correction in 
a He-like ion for the (Isf)j configuration (in case of mixing configurations, a diagonal matrix 
element should be taken), and is the total angular momentum of the v electron. By employing 
the identity (ITOb . we check that our numerical results for He-like ions are in a very good agreement 
with our previous calculations for Li-like ions lioi . 

Our calculations of the screened self-energy and two-photon exchange corrections, combined 
with the results for the screened vacuum-polarization from Ref. IbtI (with the off-diagonal matrix 
elements corrected in this paper, see below), complete the evaluation of the QED correction to 
first order inl/Z and to all orders in aZ for n = 2 states of He-like ions. As is known, the aZ 
expansion of two-electron QED effects starts with a'^{aZy. The two-photon exchange correction 
contains also contributions of previous orders in aZ that can be derived from the Breit equation. 
We separate the "pure" QED part of the two-photon exchange contribution (AE^^) as 

AE2ph = «'[ao + {aZya2] + AE^^^ , (71) 



where Ai?2ph is the total two-photon exchange correction and AE^^ contributes to order 
a'^{aZ)^ and higher. In order to extract numerical values for AE^^ from our results for Ai?2ph 
without losses in accuracy, accurate values for the coefficients ao and 02 are needed. We calculate 
them by fitting our results for the two-photon exchange correction obtained within many-body per- 
turbation theory. A large number of fitting points and inclusion of fraction values for the nuclear 
charge number (up to Z = 0.1) allowed us to achieve better accuracy than in previous calculations 
of similar coefficients (e.g., Refs. Inllsollsii '). The numerical results for the coefficients ag and 02 
for all states under consideration are tabulated in the second and in the third column of Table HiH 
respectively. 

In Table |lVl we collect all two-electron QED contributions for n = 1 and n = 2 states of 
He-like ions. The screened self-energy and two-photon exchange corrections are calculated in 
the present work; in the table they are labeled as "Scr.SE" and "2-ph.exch.", respectively. ^The 
screened vacuum-polarization correction was first evaluated in our previous investigation [37]. In 
the present work, we correct an error made in Ref. IbtI for the off-diagonal matrix element and 
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extend our calculation to the region 10 < Z < 20. Numerical values for the screened vacuum- 
polarization correction are listed in Table II VI under entry "Scr.VP". 

1 orders in aZ can be compared 



Our results for the two-electron QED correction calculated to a 



with the results obtained within the aZ expansion, which reads 



m 



asi In aZ + a^o + (aZ) G2^'{aZ) 



(72) 



Ref. [1] and numerical results for the two-electron Bethe logarithms 
expansion coefficients of expectation values of various operators IitI IsS 



where the function G2ei'{ceZ) is the higher-order remainder that is not known analytically at 
present We obtain numerical values for the coefficients 031 and 039 b;^sing formulas from 

qM and for the 1/Z- 
The only coefficient 

whose numerical value was not available in the literature was the anomalous-magnetic moment 
correction for the off-diagonal matrix element. This is the first-order 1 /Z-expansion term of the 

n 

matrix element of the operator a/n^H'^' + ifg") (see E^s. (27) and (28) of Ref. n53 ). The result 
of our calculation of this correction (denoted in Ref. 1 17] as A_Eanom) for the off-diagonal term in 
the LS coupling reads 

AEf„^^(offdiag) = a^{aZf 0.010110 . (73) 

Numerical values for the coefficients 031 and 030 for all states under consideration are listed in the 
third and in the fourth columns of Table Hill respectively. 

In Fig. m we plot our numerical results together with the contribution of the first two terms 
of the aZ expansion (dashed line). In addition, we also plot the two-electron QED contribution, 

n n 

as evaluated by Drake |17I] (dotted line). It was obtained according Eqs. (2)-(9) of Ref. II17I1 . 
keeping the contribution of first order in 1/Z only. (We note that Eq. (8) of Ref. IitI contains 
a misprint; its right-hand-side should be multiplied by Z.) Expressions obtained in this way are 
exact to the leading order a^(aZ)^. They also contain some higher-order contributions, due to 
all-order results for the one-electron QED correction employed for the evaluation of the Ei^i term 
(Eq. (2) of Ref. |17]). We observe a good agreement of our results with the previously known 
contributions and conclude that Drake's values fall much closer to our all-order results than the 
pure aZ-expansion contribution. 

For mixing states (ls2pi/2)i and (Is2p3/2)i, Fig.HJpresents a comparison for individual diago- 
nal and off-diagonal matrix elements. It should be mentioned that, generally speaking, comparison 
of different methods should be performed for the physical energies obtained after the diagonaliza- 
tion of the total matrix and not for the individual matrix elements, since matrices with the same 
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eigenvalues can differ by a unitary transformation. We see from Fig. |4l however, that our re- 
sults are in a good agreement with the aZ-expansion contributions also for the individual matrix 
elements. 

An agreement found with the leading term of the olZ expansion offers us a possibility to ob- 
tain the next-to-leading contribution, which is not known analytically at present, and in this way 
to extend the results of our calculations to lower values of Z. We thus isolate the higher-order 
remainder Ggei ' ('^^) [^^^ Eq. (l72t l from our numerical data and fit it to the form 

G^^-(aZ) = 041 Ina^ + 040 + • •) • (74) 

Fitted values for the coefficients 041 and 049 are presented in the last two columns of Table HiH It 
should be stressed that these coefficients were obtained in the jj-coupling scheme with the wave 
functions defined in case of mixing states by Eqs. ([T]), ©. 

There is a way to check the self-consistency of the numerical results for individual matrix 
elements, which allows us to check each two-electron QED contribution separately. We note that, 
in the LS coupling, the only contribution to the off-diagonal matrix element to order (y^[oLZ)^ is 
that of the anomalous magnetic moment correction Ai^anom^ Eq. W^ . Therefore, for the two- 
photon exchange and screened vacuum-polarization corrections, the off-diagonal matrix element 
in the LS coupling is zero. In this case, the following identity is valid in the jj-coupling scheme 
(to the order a^(aZ)^) 

\/2[AE(i.2p,/,). - Ai?(i.2p3/.)J = - Ai^f^diag , (75) 

where Ai?, stand for the corresponding matrix elements. For the screened self-energy correction, 
the off-diagonal matrix element in the LS coupling (Ai^^j^jj^g) is nonzero and the corresponding 
identity reads 

V^[Ai?(ls2p,/,), - AE(i.2p3/.) J + Ai^^'Liag = 3AE„^4,g . (76) 

Fulfillment of these identities for individual two-electron QED contributions is checked in TablelVl 
For the screened self-energy and vacuum-polarization correction, the fulfillment is obvious from 
the table. For the two-photon exchange correction, the difference between the right- and left-hand- 
side is very close to 3{aZ)^ eV in all cases listed and, therefore, should be ascribed to higher-order 
contributions, for which the identity is not valid anymore. 



21 



III. ENERGIES OF n = 1 AND n = 2 STATES OF HE-LIKE IONS 



In this section we collect all contributions available to the ionization energies of n = 1 and 
n = 2 states of He-like ions. Individual corrections for selected ions are listed in Table IVTl A 
description of contributions presented there is given below. 

Dirac energy. A^EDirac is the Dirac value for the ionization energy of the valence electron 
including the finite-nuclear-size effect. The energy levels were calculated employing the two- 
parameter Fermi model for the nuclear-charge distribution. Parameters of the Fermi model were 



expressed in terms of the root-mean- square (rms) radius (see, e.g., Ref. mM), whose actual values 
were taken from Refs. |55, 3.0. Si- For each value of Z, the nuclear parameters for the isotope 
with the largest abundance (with the longest life time) were chosen. An approximate formula from 
Ref. [59] was employed for calculating rms radii for ions with no experimental data available. In 
the table, we present also an estimation of the uncertainty of the nuclear-size effect. In all cases 
except Z = 80, 82, 83, 90, and 92, this uncertainty was evaluated by taking the one-percent 
variation of the rms radius. For the above mentioned exceptions, the rms-radii are supposed to 
be known more precisely. In our calculation we employed the following values: 5.467(6) Fm for 
Z = 80, 5.504(25) Fm for Z = 82, 5.533(20) Fm for Z = 83, 5.802(4) Fm for Z = 90, and 
5.860(2) Fm for Z = 92. The uncertainty of the nuclear-size effect in these cases was evaluated 
by adding quadratically two errors, one obtained by varying the rms radius within the error bars 
given and the other obtained by changing the model of the nuclear-charge distribution (the Fermi 
and the homogeneously-charged- sphere model were employed). 

Electron-electron interaction correction. AEiy^t incorporates corrections that can be derived 
from the Breit equation. It consists of 3 parts, 

which correspond to the one, two, and three and more photon exchange, respectively. In notations 



of SecUl the one-photon exchange correction is written as 11331 



ilPi2 fcifc2 

(A, 



(78) 



p 

where Ai = Spi-^ — Sk^ and A2 = epi^ — €k2 - Its numerical evaluation was carried out employing 
the Fermi model for the nuclear-charge distribution; accurate numerical results for this correction 
can be found in Ref. IstI . AEfp^'* represents the two-photon exchange correction within the 

22 



a [a 



Zy approximation and is given by the first two terms in Eq. (TtTT) . with the coefficients oq 



and a2 listed in Table Unj The contribution due to the exchange by three and more photons was 
evaluated by summing terms of the l/Z expansion, with the corresponding coefficients taken from 
Refs. H Isli for the nonrelativistic energy and from Ref. jivl] for the Breit-Pauli correction. 

One-electron QED correction. AE^^^ is the sum of the one-loop and two-loop one-electron 
QED corrections. The one-loop self-energy correction for Is, 2s, and 2pi/2 states and Z > 26 
(including the nuclear-size effect) was tabulated in Ref. [60] by using the method developed by 
Moh. and co-worke. U □ □ . For lower values of Z and for Che 2r,, s,a,e, we used aconrbi- 
nation of our own calculation and an interpolation of the point-nucleus results from Ref. |63|] . The 
Uehling part of the one-loop vacuum-polarization correction was calculated in this work for the 
Fermi nuclear model. The Wichmann-KroU part of the vacuum-polarization correction was tabu- 
lated for Z > 30 in Ref. [64]. For lower values of Z, it was calculated in this work by employing 
the asymptotic-expansion formulas for the Wichmann-KroU potential ll65h . 

The two-loop one-electron QED correction is calculated to all orders in aZ only for the Is 
state up to now, see Ref. [48] and references therein. For excited states, one has to rely on the aZ 
expansion, which reads (see review references .erern, and .o.recencscud,es0y, 



L^Bq2 + L Bqi + G 



h.o. I 
"^210 ' 



a 



(79) 



where L = \n.[{aZ) G\[°-{aZ) = Bqq + {aZ){- ■ ■ ) is the higher-order remainder, and the 



coefficients Bij are 
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BUns) = Y - ln(2n) + ^ - i + V^(n) + , (83) 
4 ri^ - 1 

B,2{np) = (84) 

(Is) = 50.344005, (85) 

(2s) = 42.447669, (86) 

56o(ls) = -61.6(9), (87) 

56o(2s) = -53.2(8), (88) 

where C, is the Riemann zeta function, ^ is the logarithmic derivative of the gamma function, and 
C = 0.577261 ... is the Euler constant. Great care should be taken employing the aZ expan- 
sion for the estimation of the total correction for middle- and high-Z ions, due to a very slow 
convergence of this expansion. In addition, it was found lately \6% that the numerical all-order 
results do not agree well with the analytical calculations to order a^(aZ)^. A possible reason for 
this disagreement lyoi can be an incompleteness of the analytical results (1551 . (15^ for the Bqi 
coefficient. 

In order to extrapolate the all-order numerical results of Ref. [ B ] to the region Z= 12-39 for 
the Is state and to estimate the two-loop correction for excited states, we separate the Is higher- 
order remainder G2[" {aZ) from the numerical data of Ref. [48]. We observe that this function 
is smoothly behaving and can be reasonably approximated by a polynomial. We thus employ a 
linear (parabolic) fit to the function Ggio' ('^^) order to extrapolate the higher-order contribution 
to the region Z=12-39. For 2s state, we employ the same values for the higher-order contribution 
and ascribe the uncertainty of 50% to them. For p states, no analytical calculations for the Bqi 
coefficient exist up to now. We thus separate from the Is numerical results of Ref. [48] the function 

{aZ)=L B,, + G\:- {aZ) , (89) 

divide it by a factor of 8, and take the result as the uncertainty for the higher-order contribution for 
p states. 

Two-electron QED correction. AE^^^ is evaluated in SecHIl the data are taken from TablellVI 
Higher-order QED correction. AE^^^ represents the contribution of QED effects of relative 
order l/Z"^ and higher. This correction was evaluated by formulas presented in Ref. [Ill sup- 
pressing terms that contribute to orders and 1/Z. Its uncertainty was obtained by taking the 
relative deviation of the QED contribution to order 1/Z calculated according to Ref. Iivll from the 
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results of its exact evaluation presented in this work. (The corresponding comparison is presented 
in Fig. El) 

Relativistic recoil correction. Ai^rec consists of the one-electron and the two-electron part. The 
one-electron relativistic recoil correction was evaluated to all orders in a Z in a series of papers 



7311 . In our compilation, we employed the finite-nucleus results of Ref. ItbI for the Is 
state, the point-nucleus results of Ref. 117 111 for the 2s and 2pi/2 states, and those of Ref. 

Q for 

the 2p3/2 state. The two-electron recoil contribution is given by the sum of the mass-polarization 
correction and the electron-electron interaction correction to the one-electron nuclear recoil. The 
nonrelativistic part of the mass-polarization correction was evaluated by summing the terms of the 

n 

l/Z expansion of the matrix element (pj - P2) taken from Ref. P17l]. The known relativistic part 
of this correction of order (aZ)^m/M ll74ll was also included. The electron-electron interaction 
correction to the one-electron nuclear recoil was taken into account in the nonrelativistic limit. It 
was estimated as (— m/M) Ai?2eh where A£^2ei is the total two-electron correction. 

In the last column of Table IvTl we present the total values for the ionization energies, which 
are given by the sum of all corrections mentioned so far. For lead, thorium, and uranium, the total 

nn 

values include also the nuclear-polarization correction 175U76I]. Analyzing the main sources of 
uncertainties listed in the table, we conclude that in the low-Z region the main error comes from 
the two-electron QED corrections, namely from the two-photon exchange contribution. In the 
high-Z region, main sources of uncertainty are the one-electron two-loop QED correction (mostly, 
the two-loop self-energy correction) and the experimental values for the rms nuclear radii. 

In Table IVIIL the total ionization energies of n = 1 and n = 2 states of He-like ions with 
Z = 12 — 100 are listed. We start our compilation with Z = 12 since this is the point where 
the new terms accounted for in our calculation (~ a^{aZY) become comparable with omitted 
higher-order effects (~ a^{aZy). 

In Fig. 121 our results are c omp ared with the theoretical values obtained previously in calcu- 
lations of different types Since our evaluation is the first one complete to the or- 
der a^{aZY, it is interesting to analyze the difference between various calculations in units of 
a^{aZY. First of all, we note a significant deviation of our values from the recent results by 
Cheng and Chen ll22ll . which arises from an incomplete treatment of QED corrections employed 
in that work. The authors evaluate the QED correction to all orders in a Z at the one-loop level, 
employing a symmetric model potential in order to account for the electron-electron interaction. 
This approximation works reasonably well in the high-Z region, but for ions with 22 < Z < 36 
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(as presented in the paper), the accuracy of this approximation turns out to be lower than that of 
Drake's approach based on the exact aZ expansion ifni . We mention that a previous investigation 
by these authors [ 19] employed the QED correction as evaluated by Drake. Its results agree well 
with those by Plante et al. ll20l] and thus are in a better agreement with our numerical values. 
For the ^Sq and 2 ^Po,i states, we observe also a distinct deviation of our ionization ener gies 



from the results by Drake A similar deviation was reported previously in the literature 1118 
where it was attributed to corrections of order a^{aZY to the electron-electron interaction 
that were not accounted for by Drake's unified method but can be (to a certain extent) included 
by methods based on the no-pair QED Hamiltonian h-regularities of the Z-dependence of 
the plotted difference, which can be observed for S states in the medium- and high-Z region, is 
explained by more recent values for the rms nuclear radii employed in the present calculation. 

As can be seen from Fig. |5l the best agreement is found with the calculation by Plante et al. 
1 20]. It is to be noted that the results by Johnson and Sapirstein lisi and by Chen et al. I19I1 



obtained by different methods but on the same level of sophistication are in a very good agreement 
with the ones by Plante and co-workers. Whereas all these results are incomplete to order {aZY, 
we conclude that the remaining contribution of this order is rather small for all n = 2 states, which 
explains a good agreement of these results with the experimental data. Only for the 1 ^Si state, we 
observe a significant new contribution of about 0.5 a^{aZY. We mention, however, that despite 
of a good agreement observed for the n = 2 states, the results by Plante et al. are well outside of 
the estimated error bars of the present theoretical values for most middle- and high-Z ions. 



In Table I Villi we list transition energies for which experimental results are available. Compari- 
son is made with the MBPT calculation by Johnson and Sapirstein [18], with the CI calculations by 
Chen et al. 13], and with the all-order many-body treatment by Plante et al. \q^ . These studies 
are, according to our analysis, the most complete ones among the previous calculations. We recall 
that in all these investigations QED corrections were taken as evaluated by Drake 0. The dif- 
ference between them, therefore, is related only to the part arising from the no-pair Hamiltonian, 
often referred to as the "structure" part. 

We observe a generally good agreement of theoretical predictions with experimental data. De- 
spite of the significant amount of available experimental information, the experimental uncertainty 
in the region of Z under consideration is generally larger than the difference between the calcula- 
tions analyzed in Table IVIlTl Among few exceptions are the recent high-precision measurements 
of the 2 ^Pi - 2 ^5*0 transition energy in silicon {Z = 14) S and the 2 ^Pq - 2 ^Pi transition 
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energy in magnesium (Z = 12) 119111 . whose accuracy is much higher than that of the theoretical 
predictions. However, at these relatively low values of Z, our treatment is basically equivalent to 
the previous studies, and the difference between the calculations can not be effectively probed in 
comparison with these measurements. When Z increases, deviation of our values from the results 
of the previous studies becomes more prominent, but the experimental uncertainty is much lower 
for higher Z. A compromise is found to be argon (Z = 18), where the experimental determination 
of the 2 ^Po,2 — 2 ^Si transition energies by Kukla et al. 11821] demonstrated a 2o deviation from 
the previous theoretical results. Our calculation brings the theoretical and experimental results 
in agreement for the 2 ^Pq — 2 ^5*1 transition and reduces the discrepancy for the 2 — 2 '^Si 
transition to 0.5 cr. 

An important feature of He-like ions is that they provide a possibility to study the effects of 
parity non-conservation lioi il|. The 2^5*0 — 2^Pq transition in He-like Eu ion {Z = 63) is 
presently considered as the best candidate for future experiments loi . The effect is enhanced by 
the fact that the 2 and 2 ^Pq levels cross each other in a vicinity of Z = 63. Another crossing 
point of the levels occurs around Z = 90 but it seems to be less promising for the experimental 
observation of the effect. In Table El we list the results of different theoretical evaluations for 
the 2^Pq — 2 transition energy in ions near the crossing points. A significant discrepancy is 
observed between different theoretical evaluations, which is due to the smallness of the energy 
difference for these ions. We mention a significant deviation of our values from the recent results 
by Andreev et al. ll4fli . In that work, the authors performed an ab initio calculation of the two- 
photon exchange correction, whose numerical values agree well with those obtained in this paper. 
However, evaluating the total transition energy, the authors used an estimation for the screened 
self-energy correction (that was not calculated at that moment), which is the main source of the 
disagreement observed. 

Summarizing, in this investigation we performed ab initio QED calculations of the screened 
self-energy correction and the two-photon exchange correction forn = 1 and n = 2 states of He- 
like ions with Z > 12. This evaluation completes the rigorous treatment of all two-electron QED 
corrections of order to all orders in aZ and significantly improves the theoretical accuracy for 
the energy values, especially in the high-Z region. Unlike all previous calculations, the results 
obtained are complete through order a^{aZY\ uncalculated terms enter through three-photon- 
exchange QED effects (~ a^{aZy' and higher) and through higher-order one-electron two-loop 
QED corrections (~ a'^{aZy and higher). 
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TABLE I: The two-photon exchange correction for n = 1 and 
n = 2 states of He-like ions, in eV. For mixing configurations, 
(ls2pi/2)i and (Is2p3/2)i stand for the diagonal matrix elements 
of the operator H [see Eqs. (fTOb l. whereas "off-diag." labels 
the off-diagonal matrix elements. 



Z (Isls)o (ls2s)o (ls2s)i (ls2pi/2)o (ls2pi/2)i {Is2py2)i (l'S2p3/2)2 off-diag. 



12 


-4.4186 


-3.1741 


-1 


.2991 


-2.0506 


-2.7789 


-3.5332 


-1.9964 


-1.0711 


14 


-4.4645 


-3.1952 


-1 


.3024 


-2.0741 


-2.7899 


-3.5413 


-2.0000 


-1.0686 






-3.19541* 


-1 


.30240* 












16 


-4.5173 


-3.2196 


-1 


.3062 


-2.1015 


-2.8027 


-3.5507 


-2.0041 


-1.0658 


18 


-4.5770 


-3.2473 


-1 


.3106 


-2.1328 


-2.8173 


-3.5613 


-2.0088 


-1.0626 






-3.24753* 


-1 


.31057* 




-2.8168'= 


-3.5603*= 




-1.0618^ 


20 


-4.6435 


-3.2784 


-1 


.3154 


-2.1682 


-2.8337 


-3.5733 


-2.0141 


-1.0589 




-4.6447« 


















28 


-4.9784 


-3.4378 


-1 


.3405 


-2.3532 


-2.9182 


-3.6340 


-2.0405 


-1.0406 


30 


-5.0795 


-3.4868 


-1 


.3483 


-2.4111 


-2.9443 


-3.6525 


-2.0484 


-1.0350 




-5.0812" 


-3.48716* 


-1 


.34827* 


-2.41112'^ 


-2.9439'= 


-3.6506^ 


-2.04834'=' 


-1.0350^ 






-3.473^ 


-1 


.348'= 


















— i, 


.o4ooo 












32 


-5.1877 


-3.5396 


-1 


.3566 


-2.4741 


-2.9725 


-3.6724 


-2.0568 


-1.0291 


40 


-5.6924 


-3.7919 


-1 


.3961 


-2.7817 


-3.1072 


-3.7658 


-2.0956 


-1.0015 




-5.6945" 




-1 


.39621'' 


-2.78172'^ 


-3.1082'= 


-3.7641^ 


-2.09545'' 


-1.0008^ 


47 


-6.2332 


-4.0719 


-1 


.4395 


-3.1351 


-3.2575 


-3.8668 


-2.1358 


-0.9724 


50 


-6.4951 


-4.2110 


-1 


.4609 


-3.3148 


-3.3323 


-3.9159 


-2.1548 


-0.9586 




-6.4975" 




-1 


.46120'^ 


-3.31489'^ 


-3.333'= 


-3.915*= 


-2.15465'=' 


-0.955'= 


54 


-6.8742 


-4.4162 


-1 


.4923 


-3.5848 


-3.4429 


-3.9871 


-2.1816 


-0.9387 


60 


-7.5114 


-4.7714 


-1 


.5459 


-4.0642 


-3.6348 


-4.1066 


-2.2251 


-0.9064 




-7.5142" 


-4.77215* 


-1 


.54587* 


-4.068'= 


-3.635*= 


-4.105'= 


-2.22510'=' 


-0.893'= 






-4.78r 


-1 


.542'= 


-4.06446'^ 











29 











-1 


54558^^ 
















66 


-8.2393 


-5 


1924 


-1 


6082 


-4.6505 


-3.8632 


-4.2426 


-2 


2724 


-0 


8708 






-5 


194"^ 


-1 


605'' 


-4.670" 














70 


-8.7812 


-5 


5159 


-1 


6552 


-5.1131 


-4.0394 


-4.3430 


-2 


3060 


-0 


8453 




-8.7847" 


-5 


515'' 


-1 


648" 


-5.117" 


-4.038'' 


-4.339" 


-2 


30573*^ 


-0 


801" 










-1 


65478'* 


-5.11403'^ 














74 


-9.3739 


-5 


8794 


-1 


7073 


-5.6441 


-4.2381 


-4.4517 


-2 


3412 


-0 


8184 


79 


-10.1957 


-6 


3996 


-1 


7803 


-6.4220 


-4.5238 


-4.5999 


-2 


3877 


-0 


7826 


80 


-10.3719 


-6 


5135 


-1 


7961 


-6.5950 


-4.5866 


-4.6312 


-2 


3974 


-0 


7752 




-10.375^* 


-6 


504^ 


-1 


789" 


-6.598" 


-4.585^ 


-4.628" 


-2 


39806'^ 


-0 


771" 










-1 


79562'^ 


-6.59593'' 














82 


-10.7375 


-6 


7524 


-1 


8289 


-6.9607 


-4.7185 


-4.6957 


-2 


4170 


-0 


7601 


83 


-10.9271 


-6 


8776 


-1 


8460 


-7.1540 


-4.7877 


-4.7288 


-2 


4270 


-0 


7524 


90 


-12.3979 


-7 


8792 


-1 


9790 


-8.7331 


-5.3458 


-4.9780 


-2 


5005 


-0 


6957 



-12.403'' 
92 -12.8714 



100 -15.0772 
-15.0805'> 



8.2122 


-2.0221 


-9 


2701 


-5.5329 


-5.0550 


-2.5228 


-0.6787 


8.21306^ 


-2.02199'' 


-9 


274" 


-5.531" 


-5.053" 


-2.52228^^ 


-0.683" 


8.184" 


-2.018" 
-2.02034'' 


-9 


27598'' 










9.8239 


-2.2223 


-11 


9330 


-6.4484 


-5.3900 


-2.6191 


-0.6058 



'^ Blundell et al. UM,^ Asen et al. MM," Andreev et al. ||39l|40|]/ Mohr and Sapirstein ||38|l, 
" Andreev et al. Q- 

TABLE IV: Two-electron QED coiTection for n = 1 and n = 2 states of He-like ions, in 
eV. 



Z State 


ScrSE 


Scr.VP 


2-ph.exch. 


Total 


Z State 


Scr.SE 


Scr.VP 


2-ph.exch. 


Total 


12 {ls)2 


-0.0405 


0.0021 


0.0031(1) 


-0.0353(1) 


60 (ls)2 


-2.4392(2) 


0.3800(1) 


0.0662(1) 


-1.9930(2) 


(ls2s)o 


-0.0088 


0.0004 


0.0004(1) 


-0.0080(1) 


(ls2s)o 


-0.6267(2) 


0.0923 


-0.1914(1) 


-0.7258(2) 


(ls2s)i 


-0.0055 


0.0003 


-0.0001(1) 


-0.0053(1) 


(ls2s)i 


-0.3377(2) 


0.0484 


-0.0327(1) 


-0.3219(2) 


(ls2pi/2)o 


-0.0017 


0.0001 


-0.0008(1) 


-0.0024(1) 


(ls2pi/2)o 


-0.1569(1) 


0.0311 


-0.4945(1) 


-0.6203(1) 


(ls2pi/2)i 


-0.0013 


0.0001 


-0.0002(1) 


-0.0014(1) 


(ls2pi/2)i 


-0.1227(1) 


0.0190 


-0.1330(1) 


-0.2367(1) 


(Is2p3/2)l 


-0.0010 


0.0000 


0.0000(1) 


-0.0010(1) 


(Is2p3/2)l 


-0.1063(2) 


0.0046 


-0.0349(1) 


-0.1366(2) 
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TABLE II: Screened self-energy correction for n = 1 and n = 2 states of He-like ions, in units of F{aZ). 
In case of mixing configurations, contributions to the matrix elements Hn. are given; labeling is as in TableU 



z 


(Isls)o 


(ls2s)o 


(ls2s)i 


(ls2pi/2)o 


(ls2pi/2)i 


(Is2p3/2)l 


(ls2p3/2)2 


off-diag. 


12 


-2.2139(8) 


-0.4841(5) 


-0.3031(5) 


-0.0917(6) 


-0.0691(6) 


-0.0556(7) 


-0.1350(7) 


0.0533(2) 


14 


-2.0543(6) 


-0.4519(4) 


-0.2821(4) 


-0.0845(5) 


-0.0646(5) 


-0.0537(6) 


-0.1266(6) 


0.0490(1) 


16 


-1.9217(3) 


-0.4248(3) 


-0.2646(3) 


-0.0783(4) 


-0.0605(4) 


-0.0517(4) 


-0.1197(4) 


0.04559(5) 


18 


-1.8097(3) 


-0.4021(3) 


-0.2496(3) 


-0.0733(2) 


-0.0571(2) 


-0.0501(2) 


-0.1137(2) 


0.04266(3) 


20 


-1.7137(3) 


-0.3828(3) 


-0.2368(3) 


-0.0693(1) 


-0.0544(1) 


-0.0488(2) 


-0.1086(2) 


0.04013(3) 


30 


— 1.3888(2) 


—0.3194(2) 


—0.1930(2) 


—0.0581(1) 


—0.0470(1) 


—0.0452(2) 


—0.0913(2) 


0.03146(2) 


40 


-1.2112(1) 


-0.2879(1) 


-0.1685(1) 


-0.05588(7) 


-0.04542(7) 


-0.0442(1) 


-0.0817(1) 


0.02639(2) 


50 


-1.1134(1) 


-0.2746(1) 


-0.1547(1) 


-0.05963(8) 


-0.04784(8) 


-0.0449(1) 


-0.0761(1) 


0.02312(2) 


60 


-1.0679(1) 


-0.2744(1) 


-0.1478(1) 


-0.06871(6) 


-0.05371(6) 


-0.0465(1) 


-0.0729(1) 


0.02087(1) 


70 


-1.06281(5) 


-0.28559(5) 


-0.14670(5) 


-0.08394(5) 


-0.06349(5) 


-0.04896(7) 


-0.07136(7) 


0.019257(8) 


80 


-1.09510(3) 


-0.30916(3) 


-0.15096(3) 


-0.10779(2) 


-0.07864(2) 


-0.05197(7) 


-0.07091(7) 


0.018047(6) 


83 


-1.11237(2) 


-0.31903(2) 


-0.15336(2) 


-0.11728(2) 


-0.08463(2) 


-0.05294(7) 


-0.07094(7) 


0.017741(5) 


90 


-1.16760(2) 


-0.34804(2) 


-0.16122(2) 


-0.14526(1) 


-0.10222(1) 


-0.05530(7) 


-0.07130(7) 


0.017104(3) 


92 


-1.18776(2) 


-0.35814(2) 


-0.16413(2) 


-0.15515(1) 


-0.10841(1) 


-0.05600(7) 


-0.07148(7) 


0.016939(3) 


100 


-1.29293(2) 


-0.40917(2) 


-0.17942(2) 


-0.20688(3) 


-0.14073(3) 


-0.05881(7) 


-0.07250(7) 


0.016343(3) 



(ls2p3/2)2 


-0.0025 


0.0001 


-0.0001(1) 


-0.0024(1) 


(ls2p3/2)2 


-0.1666(2) 


0.0159 


0.0081(1) 


-0.1426(2) 


off-diag. 


0.0010 


-0.0001 


0.0001(1) 


0.0010(1) 


off-diag. 


0.0477 


-0.0092 


0.0053(2) 


0.0437(2) 


14 (ls)2 


-0.0596 


0.0034 


0.0046(1) 


-0.0516(1) 


70 (ls)2 


-3.8548(1) 


0.7130(2) 


-0.0164(1) 


-3.1581(2) 


(ls2s)o 


-0.0131 


0.0007 


0.0005(1) 


-0.0119(1) 


(ls2s)o 


-1.0358(1) 


0.1819 


-0.4071(1) 


-1.2610(2) 


(ls2s)i 


-0.0082 


0.0005 


-0.0002(1) 


-0.0079(1) 


(ls2s)i 


-0.5321(1) 


0.0892 


-0.0615(1) 


-0.5043(2) 


(ls2pi/2)o 


-0.0025 


0.0002 


-0.0015(1) 


-0.0037(1) 


(ls2pi/2)o 


-0.3044(1) 


0.0667 


-0.9717(1) 


-1.2094(2) 


(ls2pi/2)i 


-0.0019 


0.0001 


-0.0004(1) 


-0.0021(1) 


(ls2pi/2)i 


-0.2303(1) 


0.0409 


-0.2656(1) 


-0.4550(2) 


(Is2p3/2)l 


-0.0016 


0.0001 


0.0000(1) 


-0.0015(1) 


(Is2p3/2)l 


-0.1776(2) 


0.0075 


-0.0688(1) 


-0.2388(2) 


(ls2p3/2)2 


-0.0037 


0.0002 


-0.0001(1) 


-0.0036(1) 


(ls2p3/2)2 


-0.2588(2) 


0.0266 


0.0164(1) 


-0.2158(2) 


off-diag. 


0.0014 


-0.0001 


0.0002(1) 


0.0015(1) 


off-diag. 


0.0698 


-0.0158 


0.0063(2) 


0.0603(2) 


16 (ls)2 


-0.0832 


0.0051 


0.0066(1) 


-0.0715(1) 


80 (ls)2 


-5.9289(1) 


1.2980(2) 


-0.2374(1) 


-4.8682(3) 


(ls2s)o 


-0.0184 


0.0011 


0.0005(1) 


-0.0168(1) 


(ls2s)o 


-1.6738(1) 


0.3520(1) 


-0.7946(1) 


-2.1164(2) 


(ls2s)i 


-0.0115 


0.0007 


-0.0003(1) 


-0.0110(1) 


(ls2s)i 


-0.8173(1) 


0.1615(1) 


-0.1093(1) 


-0.7652(2) 


(ls2pi/2)o 


-0.0034 


0.0003 


-0.0025(1) 


-0.0056(1) 


(ls2pi/2)o 


-0.5836(1) 


0.1429 


-1.7938(1) 


-2.2345(1) 


(ls2pi/2)i 


-0.0026 


0.0002 


-0.0006(1) 


-0.0030(1) 


(ls2pi/2)i 


-0.4258(1) 


0.0879 


-0.4988(1) 


-0.8367(1) 


(Is2p3/2)l 


-0.0022 


0.0001 


0.0000(1) 


-0.0022(1) 


(Is2p3/2)l 


-0.2814(3) 


0.0120 


-0.1232(1) 


-0.3926(3) 


(ls2p3/2)2 


-0.0052 


0.0003 


-0.0002(1) 


-0.0051(1) 


(ls2p3/2)2 


-0.3839(3) 


0.0428 


0.0279(1) 


-0.3132(3) 


off-diag. 


0.0020 


-0.0002 


0.0003(1) 


0.0021(1) 


off-diag. 


0.0977 


-0.0260 


0.0072(2) 


0.0789(2) 


18 (ls)2 


-0.1116 


0.0072 


0.0091(1) 


-0.0953(1) 


83 (ls)2 


-6.7256(1) 


1.5500(7) 


-0.3460(1) 


-5.5216(7) 


(ls2s)o 


-0.0248 


0.0015 


0.0004(1) 


-0.0228(1) 


(ls2s)o 


-1.9289(1) 


0.4286(2) 


-0.9599(1) 


-2.4602(2) 


(ls2s)i 


-0.0154 


0.0010 


-0.0004(1) 


-0.0148(1) 


(ls2s)i 


-0.9273(1) 


0.1927(2) 


-0.1289(1) 


-0.8635(2) 


(ls2pi/2)o 


-0.0045 


0.0004 


-0.0039(1) 


-0.0080(1) 


(ls2pi/2)o 


-0.7091(1) 


0.1799(1) 


-2.1377(1) 


-2.6669(1) 


(ls2pi/2)i 


-0.0035 


0.0003 


-0.0010(1) 


-0.0042(1) 


(ls2pi/2)i 


-0.5117(1) 


0.1109(1) 


-0.5977(1) 


-0.9985(1) 


(Is2p3/2)l 


-0.0031 


0.0001 


-0.0001(1) 


-0.0031(1) 


(Is2p3/2)l 


-0.3201(4) 


0.0136 


-0.1446(1) 


-0.4511(4) 
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TABLE III: Coefficients of the aZ expansion of the second-order two-electron contribution to the energy 
levels of He-like ions. In case of mixing configurations, contributions to the matrix elements are given; 



labeling is as in Table U 





iaZf 

\ } 


iaZf 


(aZf \n.aZ 


iaZf 


(qZ)^ InaZ 


(aZ)^ 

\ } 


(Isls)o 


-0.157662 


-0.6302 


1.3191 


1.6588 


0.75(15) 


-2.41(40) 


(ls2s)o 


-0.114509 


-0.2807 


0.2755 


0.3255 


0.11(2) 


-0.81(5) 


(ls2s)i 


-0.047409 


-0.0428 


0.1795 


0.1911 


0.056(11) 


-0.40(3) 


(ls2pi/2)o 


-0.072999 


-0.3035 


0.0730 


0.1063 





-0.64(2) 


(ls2pi/2)l 


-0.101008 


-0.1444 


0.0465 


0.0578 





-0.22(1) 


(Is2p3/2)l 


-0.129018 


-0.1075 


0.0201 


0.0058 





-0.10 


(ls2p3/2)2 


-0.072999 


-0.0473 


0.0730 


0.0595 


0.01(1) 


-0.14(2) 


off-diag. 


-0.039611 


0.0319 


-0.0374 


-0.0432 


-0.01(1) 


0.08(4) 



(1 ti'?T)o /nil 


—0.0070 


0.0004 


— 0002('1 ~) 


— 0068('l 1 




(1 'I'^TJo /nil 

V-'-^^f 3/2 


—0 4289('4') 


0.0489 


0318('Tl 


—0 3482('41 


off-diag. 


0.0026 


-0.0002 


0.0004(1) 


0.0028(1) 




off-diag. 


0.1073 


-0.0299 


0.0074(4) 


0.0848(4) 


20 {ls)2 


-0.1450 


0.0099 


0.0119(1) 


-0.1231(1) 


90 


{Xsf 


-9.0006(1) 


2.338(1) 


-0.7109(1) 


-7.373(1) 


(ls2s)o 


-0.0324 


0.0021 


0.0003(1) 


-0.0300(1) 




(ls2s)o 


-2.6829(1) 


0.6810(2) 


-1.4689(1) 


-3.4708(3) 


(ls2s)i 


-0.0200 


0.0014 


-0.0006(1) 


-0.0192(1) 




(ls2s)i 


-1.2428(1) 


0.2921(2) 


-0.1869(1) 


-1.1376(2) 


(ls2pi/2)o 


-0.0059 


0.0006 


-0.0059(1) 


-0.0111(1) 




(ls2pi/2)o 


-1.1197 


0.3112(2) 


-3.1842(1) 


-3.9928(2) 


(ls2pi/2)i 


-0.0046 


0.0004 


-0.0015(1) 


-0.0057(1) 




(ls2pi/2)i 


-0.7879 


0.1929(1) 


-0.9024(1) 


-1.4974(1) 


(Is2p3/2)l 


-0.0041 


0.0001 


-0.0002(1) 


-0.0042(1) 




(Is2p3/2)l 


-0.4263(5) 


0.0176 


-0.2051(1) 


-0.6138(5) 


(ls2p3/2)2 


-0.0092 


0.0006 


-0.0003(1) 


-0.0089(1) 




(ls2p3/2)2 


-0.5496(5) 


0.0663 


0.0413(1) 


-0.4420(5) 


off-diag. 


0.0034 


-0.0003 


0.0005(1) 


0.0036(1) 




off-diag. 


0.1318 


-0.0410 


0.0082(4) 


0.0991(4) 


30 


-0.3965 


0.0348 


0.0325(1) 


-0.3292(1) 


92 


{Xsf 


-9.7800(1) 


2.630(2) 


-0.8520(1) 


-8.002(2) 


(ls2s)o 


-0.0912 


0.0076 


-0.0048(1) 


-0.0884(1) 




(ls2s)o 


-2.9489(1) 


0.7770(4) 


-1.6540(1) 


-3.8259(4) 


(ls2s)i 


-0.0551 


0.0048 


-0.0024(1) 


-0.0527(1) 




(ls2s)i 


-1.3514(1) 


0.3287(2) 


-0.2074(1) 


-1.2301(3) 


(ls2pi/2)o 


-0.0166 


0.0022 


-0.0289(1) 


-0.0433(1) 




(ls2pi/2)o 


-1.2775 


0.3647(2) 


-3.5612(1) 


-4.4740(3) 


(ls2pi/2)i 


-0.0134 


0.0013 


-0.0074(1) 


-0.0195(1) 




(ls2pi/2)i 


-0.8927 


0.2262(2) 


-1.0133(1) 


-1.6798(2) 


(Is2p3/2)l 


-0.0129 


0.0005 


-0.0016(1) 


-0.0140(1) 




(Is2p3/2)l 


-0.4611(5) 


0.0188 


-0.2254(1) 


-0.6677(5) 


(ls2p3/2)2 


-0.0261 


0.0019 


-0.0003(1) 


-0.0245(1) 




(ls2p3/2)2 


-0.5886(5) 


0.0721 


0.0440(1) 


-0.4725(5) 


off-diag. 


0.0090 


-0.0010 


0.0013(1) 


0.0093(1) 




off-diag. 


0.1395 


-0.0448 


0.0084(4) 


0.1031(4) 


40 {\sf 


-0.8197 


0.0887 


0.0589(1) 


-0.6721(1) 


100 


(ls)2 


-13.6716(1) 


4.248(4) 


-1.6551(1) 


-11.079(4) 


(ls2s)o 


-0.1948 


0.0199 


-0.0252(1) 


-0.2002(1) 




(ls2s)o 


-4.3266(1) 


1.3404(8) 


-2.6409(1) 


-5.6271(8) 


(ls2s)i 


-0.1141 


0.0118 


-0.0068(1) 


-0.1091(1) 




(ls2s)i 


-1.8972(1) 


0.5366(5) 


-0.3124(1) 


-1.6730(5) 


(ls2pi/2)o 


-0.0378 


0.0060 


-0.0916(1) 


-0.1234(1) 




(ls2pi/2)o 


-2.1876(3) 


0.7067(5) 


-5.5484(1) 


-7.0293(6) 


(ls2pi/2)i 


-0.0307 


0.0036 


-0.0239(1) 


-0.0510(1) 




(ls2pi/2)i 


-1.4881(3) 


0.4408(5) 


-1.6074(1) 


-2.6547(6) 


(Is2p3/2)l 


-0.0299 


0.0012 


-0.0058(1) 


-0.0345(1) 




(Is2p3/2)l 


-0.6219(7) 


0.0234 


-0.3211(1) 


-0.9195(7) 


(ls2p3/2)2 


-0.0553 


0.0044 


0.0005(1) 


-0.0503(1) 




(ls2p3/2)2 


-0.7666(7) 


0.1009(1) 


0.0529(1) 


-0.6128(7) 


off-diag. 


0.0179 


-0.0025 


0.0025(1) 


0.0179(1) 




off-diag. 


0.1728 


-0.0630 


0.0104(4) 


0.1202(4) 
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TABLE V: Right-hand-side (r.h.s.) and left-hand-side (l.h.s.) of Eq. (I76t (for the screened self-energy cor- 
rection) and those of Eq. ( 1751 (for the screened vacuum-polarization and two-photon exchange corrections), 
in eV. The comparison is valid to the leading order in aZ only. The last column demonstrates that the 
difference (r.h.s.-l.h.s.) for the two-photon exchange correction arises predominantly from effects to order 



z 


Scr.SE 


Scr.VP 




2-ph.exch 






l.h.s. 


r.h.s. 


l.h.s. 


r.h.s. 


l.h.s. 


r.h.s. 


(r.h.s.-l.h.s.)/(aZ)4 


12 


0.0006 


0.0006 


0.0001 


0.0001 


-0.0003 


-0.0001 


3. 


14 


0.0010 


0.0009 


0.0001 


0.0001 


-0.0005 


-0.0002 


3. 


16 


0.0014 


0.0013 


0.0002 


0.0002 


-0.0008 


-0.0003 


3.0 


18 


0.0020 


0.0019 


0.0002 


0.0002 


-0.0012 


-0.0004 


3.0 


20 


0.0027 


0.0026 


0.0003 


0.0003 


-0.0018 


-0.0005 


3.0 


30 


0.0083 


0.0087 


0.0012 


0.0010 


-0.0083 


-0.0013 


3.0 


40 


0.0167 


0.0205 


0.0035 


0.0025 


-0.0256 


-0.0025 


3.18 





-1.4717(1) 


0.1920 


0.0781(1) 


-1.2016(1) 


{ls2s)o 


-0.3630(1) 


0.0446 


-0.0783(1) 


-0.3966(1) 


(ls2s)i 


-0.2044(1) 


0.0250 


-0.0159(1) 


-0.1953(1) 


(ls2pi/2)o 


-0.0788(1) 


0.0141 


-0.2289(1) 


-0.2936(1) 


{ls2pi/2)i 


-0.0632(1) 


0.0086 


-0.0606(1) 


-0.1152(1) 


{Is2p3/2)l 


-0.0593(1) 


0.0025 


-0.0156(1) 


-0.0725(1) 


{ls2p3/2)2 


-0.1006(1) 


0.0088 


0.0031(1) 


-0.0887(1) 


off-diag. 


0.0306 


-0.0050 


0.0039(1) 


0.0294(1) 



TABLE VI: Individual contributions to the ionization energies of He-like ions (with the 
opposite sign), in eV. For mixing configurations, contributions to the matrix elements are 
listed. 



Z State 






Icl 




n.o. 




Total 


12 {\sf 


-1962.9887 


200.8973 


0.2801 


-0.0353(1) 


0.0008 


0.0412 


-1761.8045(1) 


(ls2s)o 


-490.9832 


72.9751 


0.0371 


-0.0080(1) 


0.0005 


0.0096 


-417.9688(1) 


(ls2s)i 


-490.9832 


60.2485 


0.0371 


-0.0053(1) 


0.0001 


0.0099 


-430.6928(1) 


(ls2pi/2)o 


-490.9834 


72.1736 


-0.0010 


-0.0024(1) 


0.0002 


0.0037 


-418.8092(1) 


(ls2pi/2)l 


-490.9834 


74.9743 


-0.0010 


-0.0014(1) 


0.0002 


0.0075 


-416.0038(1) 


(Is2p3/2)l 


-490.0399 


77.7699 


0.0012 


-0.0010(1) 


0.0001 


0.0113 


-412.2584(1) 


(ls2p3/2)2 


-490.0399 


71.7742 


0.0012 


-0.0024(1) 


0.0003 


0.0036 


-418.2630(1) 


off-diag. 





4.1676 





0.0010(1) 


-0.0001 


0.0054 


4.1739(1) 
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14 




-2673.7078 


235.5743 


0.4778 


-0.0516(1) 


0.0008 


0.0487 


-2437.6577(1) 




(ls2s)o 


-668.8650 


85.8200 


0.0637 


-0.0119(1) 


0.0006 


0.0115 


-582.9812(1) 




(ls2s)i 


-668.8650 


70.5889 


0.0637 


-0.0079(1) 


0.0002 


0.0118 


-598.2083(1) 




(ls2pi/2)o 


-668.8654 


84.7685 


-0.0017 


-0.0037(1) 


0.0002(1) 


0.0044 


-584.0977(1) 




(ls2pi/2)i 


-668.8654 


88.0659 


-0.0017 


-0.0021(1) 


0.0002 


0.0090 


-580.7941(1) 




(Is2p3/2)l 


-667.1144 


91.3604 


0.0022 


-0.0015(1) 


0.0001 


0.0135 


-575.7397(1) 




(ls2p3/2)2 


-667.1144 


84.1211 


0.0022 


-0.0036(1) 


0.0004 


0.0043 


-582.9899(1) 




off-diag. 





5.0013 





0.0015(1) 


-0.0002 


0.0065 


5.0091(1) 


16 




-3495.004,3 


270.4823 


0.7562 


-0.0715(1) 


0.0009 


0.0563 


-3223.7801(1) 




(ls2s)o 


-874.5000 


98.7466 


0.1014 


-0.0168(1) 


0.0007 


0.0134 


-775.6547(1) 




(ls2s)i 


-874.5000 


80.9665 


0.1014 


-0.0110(1) 


0.0002 


0.0137 


-793.4291(1) 




(ls2pi/2)o 


-874.5006 


97.4677 


-0.0028 


-0.0056(1) 


0.0003(2) 


0.0051 


-777.0359(2) 




(ls2pi/2)i 


-874.5006 


101.2216 


-0.0028 


-0.0030(1) 


0.0003 


0.0105 


-773.2742(1) 




(l«2p3/2)l 


-871.5074 


104.9766 


0.0038 


-0.0022(1) 


0.0001 


0.0158 


-766.5134(1) 




(ls2p3/2)2 


-871.5074 


96.4857 


0.0038 


-0.0051(1) 


0.0005 


0.0051 


-775.0175(1) 




off-diag. 





5.8246 





0.0021(1) 


-0.0002 


0.0076 


5.8341(1) 


18 




-4427.4152(1) 


305.6561 


1.1310(1) 


-0.0953(1) 


0.0009 


0.0575 


-4120.6651(2) 




(ls2s)o 


-1108.0563 


111.7675 


0.1525 


-0.0228(1) 


0.0007(1) 


0.0138 


-996.1445(1) 




(ls2s)i 


-1108.0563 


91.3873 


0.1525 


-0.0148(1) 


0.0003 


0.0141 


-1016.5169(1) 




(ls2pi/2)o 


-1108.0574 


110.2884 


-0.0043 


-0.0080(1) 


0.0003(3) 


0.0053 


-997.7757(3) 




(ls2pi/2)l 


-1108.0574 


114.4515 


-0.0043 


-0.0042(1) 


0.0003 


0.0108 


-993.6035(1) 




(Is2p3/2)l 


-1103.2520 


118.6221 


0.0062 


-0.0031(1) 


0.0001 


0.0162 


-984.6105(1) 




(ls2p3/2)2 


-1103.2520 


108.8712 


0.0062 


-0.0068(1) 


0.0005 


0.0052 


-994.3756(1) 




off-diag. 





6.6353 





0.0028(1) 


-0.0002 


0.0078 


6.6456(1) 


20 


(Isf 


-5471.5558(2) 


341.1317 


1.6179(2) 


-0.1231(1) 


0.0008 


0.0715 


-5128.8570(3) 




(ls2s)o 


-1369.7265 


124.8955 


0.2195(1) 


-0.0300(1) 


0.0008(1) 


0.0172 


-1244.6235(2) 




(ls2s)i 


-1369.7265 


101.8571 


0.2195(1) 


-0.0192(1) 


0.0003 


0.0175 


-1267.6513(1) 




(ls2pi/2)o 


-1369.7284 


123.2478 


-0.0063 


-0.0111(1) 


0.0004(4) 


0.0066 


-1246.4910(4) 




(ls2pi/2)l 


-1369.7284 


127.7658 


-0.0063 


-0.0057(1) 


0.0003(1) 


0.0135 


-1241.9608(1) 




(Is2p3/2)l 


-1362.3853 


132.3004 


0.0097 


-0.0042(1) 


0.0001 


0.0203 


-1230.0590(1) 




(ls2p3/2)2 


-1362.3853 


121.2809 


0.0097 


-0.0089(1) 


0.0006 


0.0065 


-1241.0965(1) 




off-diag. 





7.4312 





0.0036(1) 


-0.0003 


0.0097 


7.4442(1) 


30 


(ls)2 


-12395.3519(21) 


524.3345 


6.305(2) 


-0.3292(1) 


-0.0002 


0.1036 


-11864.9380(26) 




(ls2s)o 


-3108..3049(2) 


192.6165 


0.883(1) 


-0.0884(1) 


0.0010(4) 


0.0254 


-2914.8679(15) 




(ls2s)i 


-3108.3049(2) 


155.1570 


0.883(1) 


-0.0527(1) 


0.0005 


0.0257 


-2952.2919(14) 




(ls2pi/2)0 


-3108.3193 


190.7532 


-0.0229(3) 


-0.0433(1) 


0.0006(15) 


0.0099 


-2917.6217(15) 




(ls2pi/2)l 


-3108.3193 


195.9801 


-0.0229(3) 


-0.0195(1) 


0.0006(5) 


0.0199 


—2912.3612(6) 




(Is2p3/2)l 


-3070.5057 


201.3174 


0.0546(3) 


-0.0140(1) 


0.0001 


0.0295 


-2869.1181(4) 




(ls2p3/2)2 


-3070.5057 


183.7941 


0.0546(3) 


-0.0245(1) 


0.0012 


0.0096 


-2886.6708(4) 




off-diag. 





11.1222 





0.0093(1) 


-0.0006 


0.0140 


11.1450(1) 


40 


(15)2 


-22253.1573(98) 


720.9169 


16.315(5) 


-0.6721(1) 


-0.0030(12) 


0.1354 


-21516.466(11) 




(ls2s)o 


-5593.9685(13) 


265.1796 


2.365(7) 


-0.2002(1) 


0.0007(5) 


0.0334 


-5326.5902(69) 




(ls2s)i 


-5593.9685(13) 


210.6536 


2.365(7) 


-0.1091(1) 


0.0008 


0.0338 


-5381.0246(69) 




(ls2pi/2)o 


-5594.0369 


264.5843 


-0.040(2) 


-0.1234(1) 


0.0008(37) 


0.0134 


-5329.6015(40) 




(ls2pi/2)l 


-5594.0369 


268.0222 


-0.040(2) 


-0.0510(1) 


0.0007(12) 


0.0260 


-5326.0786(20) 




{Is2p3/2)l 


-5471.5704 


271.7376 


0.193(2) 


-0.0345(1) 


0.0001 


0.0378 


-5199.6369(16) 
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(ls2p3/2)2 


-5471.5704 


247.3427 


0.193(2) 


-0.0503(1) 


0.0018 


0.0125 


-5224.0712(16) 


off-diag. 





14.1717 





0.0179(1) 


-0.0010(1) 


0.0178 


14.2065(1) 


50 


-35226.611(37) 


936.5564 


33.961(8) 


-1.2016(1) 


-0.0077(50) 


0.1659 


-34257.137(38) 


(ls2s)o 


-8884.0997(51) 


344.7875 


5.118(22) 


-0.3966(1) 


-0.0001 


0.0412 


-8534.550(22) 


(ls2s)i 


-8884.0997(51) 


269.3108 


5.118(22) 


-0.1953(1) 


0.0011(1) 


0.0415 


-8609.824(22) 


(ls2pi/2)o 


-8884.3678(1) 


347.6605 


-0.006(5) 


-0.2936(1) 


0.0009(77) 


0.0170 


-8536.9893(91) 


(ls2pi/2)i 


-8884.3678(1) 


345.6503 


-0.006(5) 


-0.1152(1) 


0.0009(26) 


0.0315 


-8538.8067(55) 


(Is2p3/2)l 


-8575.5139 


344.0793 


0.523(5) 


-0.0725(1) 


-0.0002(2) 


0.0447 


-8230.9398(49) 


(ls2p3/2)2 


-8575.5139 


312.2766 


0.523(5) 


-0.0887(1) 


0.0025(1) 


0.0152 


-8262.7855(49) 


off-diag. 





16.3734 





0.0294(1) 


-0.0015(3) 


0.0206 


16.4220(3) 


60 


-51577.89(11) 


1178.1908 


61.92(2) 


-1.9930(2) 


-0.014(15) 


0.2152 


-50.3,39.58(12) 


(ls2s)o 


-13062.076(17) 


434.2620 


9.74(5) 


-0.7258(2) 


-0.0014(24) 


0.0537 


-12618.746(55) 


(ls2s)i 


-13062.076(17) 


332.3354 


9.74(5) 


-0.3219(2) 


0.0015(3) 


0.0541 


-12720.265(55) 


(ls2pi/2)o 


-13062.9663(6) 


443.7641 


0.20(1) 


-0.6203(1) 


0.001(14) 


0.0227 


-12619.596(19) 


(ls2pi/2)l 


-13062.9663(6) 


431.1435 


0.20(1) 


-0.2367(1) 


0.0010(46) 


0.0400 


-12631.816(13) 


(Is2p3/2)l 


-12395.4629 


418.8977 


1.20(1) 


-0.1.366(2) 


-0.0005(9) 


0.0546 


-11975.449(12) 


(ls2p3/2)2 


-12395.4629 


378.9465 


1.20(1) 


-0.1426(2) 


0.0034(1) 


0.0193 


-12015.437(12) 


oif-diag. 





17.5341 





0.0437(2) 


-0.0021(5) 


0.0246 


17.6003(6) 


70 


-71678.25(34) 


1454.7177 


103.45(5) 


-3.1581(2) 


-0.023(36) 


0.2612 


-70123.00(34) 


(ls2s)o 


-18247.262(53) 


537.4108 


17.07(10) 


-1.2610(2) 


-0.0033(84) 


0.0658 


-17693.98(11) 


(ls2s)i 


-18247.262(53) 


401.3107 


17.07(10) 


-0.5043(2) 


0.0019(5) 


0.0662 


-17829.31(11) 


(ls2pi/2)o 


-18250.1817(30) 


558.0871 


0.84(3) 


-1.2094(2) 


0.001(25) 


0.0283 


-17692.431(37) 


(ls2pi/2)l 


-18250.1817(30) 


527.6354 


0.84(3) 


-0.4550(2) 


0.0011(74) 


0.0469 


-17722.109(28) 


(Is2p3/2)l 


-16948.0254 


496.7905 


2.44(3) 


-0.2388(2) 


-0.0010(24) 


0.0608 


-16448.971(27) 
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-16948.0254 


447.7121 


2.44(3) 


-0.2158(2) 


0.0043(1) 


0.0225 
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off-diag. 





17.4750 





0.0603(2) 


-0.0030(9) 


0.0266 


17.5590(9) 


80 


-96061.17(14) 


1778.3460 


162.76(12) 


-4.8682(3) 


-0.035(78) 


0.3326 


-94124.63(20) 


(ls2s)o 


-24612.823(24) 


659.7043 


28.31(16) 


-2.1164(2) 


-0.006(21) 


0.0853 


-23926.84(16) 
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-24612.823(24) 


478.4429 


28.31(16) 
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0.0858 


-24106.74(16) 
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-24621.4092(19) 
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2.41(5) 


-2.2345(1) 


0.001(41) 


0.0369 
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639.7159 
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4.56(5) 


-0.3926(3) 


-0.0017(52) 


0.0686 
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-22253.6733 


518.9620 


4.56(5) 
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0.0289 
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32.69(18) 
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-26799.9095(69) 
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5.42(6) 
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0.0057 
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0.0297 


15.412(2) 
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245.28(27) 
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809.3517 


45.19(21) 


-3.4708(3) 


-0.009(45) 


0.1168 


-31562.77(23) 
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45.19(21) 


-1.1376(2) 
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(ls2pi/2)o 


-32440.5502(70) 
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5.82(9) 


-3.9928(2) 


0.002(66) 
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7.93(9) 
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{ls2p3/2)2 


-28337.2409 


593.1280 


7.93(9) 


-0.4420(5) 


0.0066 


0.0314 


-27736.589(93) 
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13.0506 





0.0991(4) 


-0.0052(23) 


0.0305 


13.173(2) 


92 


-132081.13(40) 


2253.9270 


265.16(33) 


-8.002(2) 


-0.05(18) 


0.4600 


-129569.84(55) 
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-34177.718(76) 


843.6057 


49.44(22) 


-3.8259(4) 


-0.009(51) 


0.1260 


-33288.42(24) 
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-34177.718(76) 


586.3549 


49.44(22) 


-1.2301(3) 


0.0030(16) 


0.1266 


-33543.06(23) 


(ls2pi/2)o 


-34211.0649(86) 


917.4965 


6.86(10) 


-4.4740(3) 


0.002(73) 


0.0531 


-33291.13(13) 


(ls2pi/2)l 


-34211.0649(86) 
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6.86(10) 


-1.6798(2) 


0.001(17) 


0.0743 
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(Is2p3/2)l 


-29649.8340 


682.1947 


8.80(10) 
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-0.003(11) 


0.0774 


-28959.44(10) 


(ls2p3/2)2 


-29649.8340 


608.3558 


8.80(10) 


-0.4725(5) 


0.0068 
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12.2592 





0.1031(4) 


-0.0054(25) 


0.0308 


12.383(3) 
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-161165.5(6.0) 


2646.5635 


358.30(63) 


-11.079(4) 


-0.06(30) 
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(ls2s)o 


-42048.7(1.2) 
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-5.6271(8) 


-0.012(86) 


0.1895 


-40984.1(1.3) 


(ls2s)i 


-42048.7(1.2) 


671.7243 


70.19(20) 


-1.6730(5) 


0.0035(23) 


0.1902 


-41308.3(1.3) 


(ls2pi/2)o 


-42127.25(19) 


1111.1289 


12.82(16) 


-7.0293(6) 


0.00(11) 


0.0759 


-41010.25(27) 


(ls2pi/2)i 


-42127.25(19) 


944.3801 


12.82(16) 


-2.6547(6) 


0.001(21) 


0.0984 


-41172.60(25) 


(Is2p3/2)l 


-35228.5685 


755.4926 


13.05(16) 


-0.9195(7) 


-0.004(16) 


0.0866 


-34460.87(16) 


(ls2p3/2)2 


-35228.5685 


670.7584 


13.05(16) 


-0.6128(7) 


0.0079(2) 


0.0382 


-34545.33(16) 


off-diag. 
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TABLE VH: Total ionization energies (in eV) for n = 1 and n = 2 states of He-like ions. 
"RMS" denotes the root-mean-square radii expressed in Fermi. 
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Po 


23 
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Pi 
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P2 
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13 
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0264(1) 
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14 
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598 
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8517(1) 


692 
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5896(1) 


16 
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17 
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3431(2) 


882 
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4785(1) 


883 
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19 
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21 
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9027(5) 


1379 
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4021(5) 


1380 
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22 
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23 
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7481 
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2019 
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3290(2) 
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26 


3.738 


8828 
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2160 


1629(7) 
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5742(7) 
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6259(10) 
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6082(4) 
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7522(2) 
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8529(2) 


27 
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1817(15) 
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8845(18) 
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8843(11) 
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5046(12) 
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2797(3) 
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3951(3) 
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3.928 


11864 


9380(26) 


2914 


8679(15) 
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2919(14) 


2917 
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2915 


0645(6) 


2866 


4147(4) 


2886 


6708(4) 


31 
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12696 


5555(31) 


3121 


9372(18) 


3160 


9268(17) 


3124 


7480(17) 


3122 


1099(7) 


3068 


3488(5) 


3089 


1461(5) 


32 


4.072 


13557 


4188(37) 


3336 


4503(21) 


3377 


0324(20) 


3339 


3118(19) 


3336 


6195(8) 


3277 


2016(5) 
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FIG. 1 : The diagram of the one-photon exchange. 




FIG. 2: The diagrams of the two-photon exchange. 
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